Evidence of Vortices on the Insulating Side of the Superconductor-Insulator Transition 
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The magnetoresistance of ultrathin insulating films of Bi has been studied with magnetic fields 
applied parallel and perpendicular to the plane of the sample. Deep in the strongly localized regime, 
the magnetoresistance is negative and independent of field orientation. As film thicknesses increase, 
the magnetoresistance becomes positive, and a difference between values measured in perpendicular 
and parallel fields appears, which is a linear function of the magnetic field and is positive. This is 
not consistent with the quantum interference picture. We suggest that it is due to vortices present 
on the insulating side of the superconductor-insulator transition. 
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In the limit of zero temperature, ultrathin films of 
metals become either insulating or superconducting de- 
pending on the strength of the disorder or the ap- 
plied magnetic field [^j. The superconductor- insulator 
(SI) transition has been described by a boson Hubbard 
Q model which also predicts a metallic phase for films 
on the margin between insulating and superconducting 
behavior whose resistance is "universal." Within this 
model, at T = 0, the superconducting state is a Cooper 
pair condensate with localized vortices, and the insulat- 
ing state is a vortex condensate with localized Cooper 
pairs. There has been no full generalization to nonzero 
temperature. The boson Hubbard model has been chal- 
lenged by the results of tunneling experiments jij carried 
out on thin films of metals which have been interpreted 
as showing that the amplitude of the order parameter is 
extremely small or zero at the SI transition and in the in- 
sulating state. On the other hand, Hall and longitudinal 
resistance studies on indium oxide films |^ , driven out of 
the superconducting state by magnetic fields, were inter- 
preted as evidence of two different insulating phases, one 
presumably a Bose insulator of localized Cooper pairs 
and the other the usual Fermi insulator of localized elec- 
trons. If this interpretation were correct, films insulating 
by virtue of disorder might behave similarly. Further- 
more, if there were vortices and Cooper pairs on the in- 
sulating side of the transition, then one might expect to 
see some evidence of vortex motion in the magnetoresis- 
tance of such films. 

The magnetoresistance (MR) of disordered two dimen- 
sional electronic systems Q in both strongly (SL) and 
weakly (WL) localized regimes is found to be dominated 
by two basic mechanisms: an orbital effect [|^,^ which is 
due to quantum interference between electron paths, and 
a spin effect |^,|l0| which arises from Zeeman splitting of 
the electron spin states. The former is flux driven, and 
therefore highly sensitive to whether the magnetic field is 
applied in the direction perpendicular to the plane of the 
sample or parallel to it, while the latter is isotropic and 
depends only on the magnetic field strength. In systems 
with strong spin-orbit (SO) coupling, both mechanisms 
predict a negative MR in the SL regime and a positive 
MR in the WL regime. The characteristic anisotropy 
due to the orbital effect, as well as the signature of the 
spin effect, has been observed in a number of experi- 
ments 1^,^. Here we report the results of MR mea- 
surements carried out on ultrathin insulating films of Bi 
close to the SI transition in magnetic fields both parallel 
and perpendicular to the plane of the sample. The ob- 
served anisotropic component of the MR is a linear func- 
tion of the magnetic field, always positive and it increases 
monotonically through the SL-WL crossover. This is not 
consistent with it being an orbital effect. We suggest that 
this contribution arises from the flow of vortices, which 
may offer new evidence in favor of the boson Hubbard 
model. 



The ultrathin Bi films used in this study were evap- 
orated on top of a lOA thick layer of amorphous Ge, 
which was pre-deposited onto the substrate which was a 
0.75mm thick single-crystal of SrTiOs (100). The sub- 
strate temperature was kept below 20 K during all depo- 
sitions and all the films were grown in situ under UHV 
conditions (~ 10~"'^°Torr). Under such circumstances, 
successive depositions could be carried out without con- 
tamination to increase the film thickness in increments of 
0.25A up to I5A. Film thicknesses were determined us- 
ing a previously calibrated quartz crystal monitor. Films 
prepared in this manner are believed to be homogeneous, 
since it has been found that they become connected at an 
average thickness on the order of one monolayer fl^ . Re- 
sistance measurements were carried out using a standard 
dc four-probe technique with current bias up to 50nA. 
Current-voltage characteristics were linear up to at least 
IfiA. Magnetic fields up to 20kG parallel to the plane of 
the sample, and up to 12kG perpendicular to the plane 
of the sample were applied using a superconducting split- 
coil magnet. 

The temperature dependence of the conductance is 
logarithmic at higher temperatures, as expected for a 
weakly localized 2D system, while at the lowest tem- 
peratures it becomes activated and could be fit by G = 
Goexp{—[TQ/T]^) where x = 1/2. At intermediate tem- 
peratures, there is a crossover region between the loga- 
rithmic and exponential behavior, which moves towards 
lower temperatures in thicker films. 

Fig. shows the MR as a function of magnetic field 
at 0.7K for the four most insulating films in both paral- 
lel and perpendicular field. The thinnest film, sample 1, 
shows a negative MR with a fractional change of up to 
12% in a 20kG field, which seems to be field-orientation 
independent. Although not shown on the graph, addi- 
tional data taken at different temperatures confirm this 
isotropic behavior. Sample 2, which is 0.25A thicker, also 
shows a negative MR, but the fractional change is now 
less than 4% in a 20kG field. For this sample, the parallel 
MR is greater than the perpendicular. This anisotropy 
becomes more obvious for sample 3, which shows a small 
positive MR for both directions of the field. The MR of 
sample 4 is even more positive and the difference between 
the parallel and the perpendicular field increases further. 

The MR changes sign from negative to positive in both 
parallel and perpendicular magnetic field at the temper- 
atures and the thicknesses for which the system is at the 
crossover between the SL and the WL regime (crossover 
being loosely defined as the region where the behavior 
of the conductance changes from activated to diffusive). 
A similar sign change has been observed in thin films of 
An, Ag and Cu by Hsu et al. [Q, where the MR was 
considered to be dominated by the orbital effect. 

In our case, the MR is dominated by a spin effect. 
One can conclude that from the fact that plotting the 
longitudinal MR against iisB/ksT makes all the data 
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collapse on a single curve. This also implies that the 
orbital contribution is negligible in parallel field. In the 
variable range hopping regime, a mechanism proposed by 
Eto ^ predicts an isotropic, negative MR due to Zeeman 
effect even in the presence of SO interactions. A model 
by Maekawa and Fukuyama also based on Zeeman 
effect, predicts a positive MR in the WL regime. Our 
results are in qualitative agreement with these pictures. 

We define the anisotropic component of the MR as 
the difference between the transverse and the longitudi- 
nal magnetoresistances Ra = (Ai?j_ — Ai?||)/i?(0) [ p^ . 
Fig. 1^ shows that Ra is a linear function of the magnetic 
field for ten different thicknesses. We will argue that this 
anisotropy cannot be due to any standard orbital effect. 
In the SL regime, in which conduction occurs through 
variable range hopping, the orbital MR can be described 
by a model originally proposed by Nguyen, Spivak and 
Shklovskii (NSS) in which an applied magnetic field 
increases the conductance. The critical percolation path 
approach |^ yields a fractional change in the resistance 
Ai?/i?(0) oc at low fields, which increases with de- 
creasing temperature. A similar behavior is predicted in 
the presence of SO scattering, although the magnitude 
of the MR is expected to decrease as the SO scattering 
increases [T^Jl7[ |. 

In a 2D system, the sample thickness is smaller than 
the hopping length, which forces all hops to be in the 
plane of the sample. This restriction on hop orientations 
leads to a MR anisotropy and a MR in a perpendicular 
field which should be much larger than in a parallel field. 
Even though the MR of sample 2 is slightly anisotropic, 
this anisotropy cannot be due to the orbital mechanism 
described above as the magnitude of the MR is smaller 
in perpendicular than in parallel fields. In addition, a 
predicted quadratic behavior with magnetic field was not 
observed. 

In the WL regime, the orbital mechanism is different. 
In systems with strong SO coupling there is a positive 
MR, which becomes larger with decreasing temperature 
as the electron phase coherence time increases. Quanti- 
tatively, the theory predicts AR/R{0) oc in the low 
field limit and AR/R{0) oc InB in the high field limit 
[^,|8|. In a 2D system, a field parallel to the plane will 
make no contribution to the fiux through the time re- 
versed loops, and no parallel MR is expected. Samples 3 
and 4 considered in Fig. Q are in the WL regime where 
the conductance shows a logarithmic dependence on tem- 
perature. The MR is positive and its magnitude is larger 
in perpendicular field, which might be expected from the 
model described above, but further analysis reveals some 
inconsistencies. Namely, it was not possible to obtain 
a satisfactory fit to the WL orbital effect theory p8| , 
which predicts Ra oc B^ in the low field. Instead, Ra is 
clearly a linear function of the magnetic field for a wide 
range of thicknesses and temperatures. Furthermore, the 
anisotropy due to the orbital effect is expected to change 



sign [|4j at the SL-WL crossover. The anisotropic re- 
sponse first occurs at the temperatures and thicknesses 
at which the films are in the SL regime, and it persists 
into the WL regime, but its magnitude is always posi- 
tive and it changes monotonically through the crossover, 
as shown on Fig. ^. All this suggests that standard or- 
bital effects cannot be the mechanism responsible for the 
anisotropy. 

Orbital MR has been observed in a number of experi- 
ments on several other systems and is often used 
to obtain various scattering times. One might therefore 
ask why we do not see any orbital effect? It is known 
that orbital MR is suppressed in the presence of strong 
SO coupling [|l6|,0, which might explain why the con- 
tribution to the MR from the orbital effect is so small 
compared to indium-oxide films. It is possible that this 
contribution would become overwhelmed by the effects 
of superconducting ordering before becoming significant 
(this would not be the case for Au, Ag and Cu films, 
which can get deep into the WL regime where the orbital 
effect is substantial, without becoming superconducting). 
Indeed, closer to the SI transition, Ra starts to deviate 
slightly from the linear field dependence (the very top of 
the Fig. |2|), acquiring a small quadratic term. This might 
be due to the orbital effect, or possibly Maki-Thompson 
(MT) effect, which has the same field dependence as the 
WL orbital effect 

If the boson-Hubbard model actually described the in- 
sulating state near the SI transition, then the insulator 
might be able to sustain point-like vortices because of the 
nonvanishing Cooper pair density. If these vortices were 



to move freely 20 1 in response to currents, they might 



produce a flux flow contribution to the MR of the sys- 
tem. This would result in Rj^-R^^ always being positive 
and proportional to the magnetic field which is the 
observed behavior. Increasing the thickness (decreasing 
the disorder) and lowering the temperature drives the 
system deeper into this putative Bose insulator state, 
(the most insulating films do not exhibit this effect and 
are Fermi insulators) where Cooper pairs and vortices 
are more likely to form. The anisotropy due to fiux fiow 
in this regime would therefore become more pronounced 
with increasing thickness and decreasing temperature, as 
observed in the measurements. In other words, the lin- 
ear response, except for its relatively small magnitude, 
very much resembles the MR due to flux flow in super- 
conducting films, which our samples become when made 
just slightly thicker. If the process were a dynamical ef- 
fect, i.e., the vortex-like MR response was occurring in 
the presence of some kind of order parameter amplitude 
and phase fluctuations, this picture might not be incom- 
patible with the interpretation of tunneling spectra on 
the insulating side of the transition, that the supercon- 
ducting gap and the average order parameter amplitude 
are both zero jij. Alternatively, and consistent with the 
original boson Hubbard picture, a nonvanishing super- 
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conducting order parameter amplitude could persist into 
the insulating regime, even though superconductivity is 
destroyed by phase fluctuations. Phase fluctuations can 
reduce the gap , as they are pair breaking 1 23 1 . This 
can lead to gapless superconductivity. It is also possi- 
ble that the local gap is destroyed by disorder, but the 
fermionic degrees of freedom are highly suppressed and 
maybe even dynamically irrelevant |20| . In support of 
this view are indications of superconducting effects in 
films of In20z driven into the insulating state by mag- 
netic fields, in the work of Paalanen et al. and in in- 
sulating films of In203 in the work of Gantmakher and 
co-workers [Q. 

Further experimental work and a theoretical model 
suitable for the intermediate disorder at finite temper- 
atures are needed to resolve this issue. In particular, we 
are attempting to detect the presence of vortices in the 
insulating regime by searching for vortex shot noise in a 
manner similar to the study of Knoedler and Voss ||2^ ] 
carried out on granular aluminum films in the supercon- 
ducting state. 
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FIG. 1. Magnetoresistance as a function of magnetic field 
for: a) four films of different nominal thicknesses: 10.50A 
(circles), 10.75A (diamonds), ll.OOA (triang les) and 11.25A 
(squares); b) sample 3 (d=11.00A) at two different tempera- 
ture: 0.7K (circles) and 0.3K (triangles). Full symbols repre- 
sent the perpendicular field, and open symbols represent the 
parallel field. 



FIG. 2. The difference between the perpendicular and par- 
allel magnetoresistance Ra as a function of magnetic field for 
ten films of different thicknesses at 0.7K. Full lines are linear 
fits. The zeros of the vertical scale are offset for clarity. 



FIG. 3. The slope of Ra vs B, which can be thought of as a 
measure of the anisotropy, as a function of film thickness for 
six different films at T=0.7K. Inset: The slope of Ra vs B as 
a function of temperature for sample 3 (d=11.00A). The lines 
are guides to the eye. The arrows show where the system is 
at the SL-WL crossover. 
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